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A pot-pourri of results in QCD from large lattice simulations on the CM5* 

Tanmoy Bhattacharya and Rajan Gupta ^ 

''T-S Group, MS B285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 U. S. A. 

We present a status report on simulations being done on 32^ X 64 lattices at /3 = 6.0 using quenclied Wilson 
fermions. Results for the spectrum, decay constants, the kaon B-parameter Bk. and semi-leptonic form factors 
are given based on the current statistical sample of 28 configurations. These "Grand Challenge" calculations are 
being done on the CM5 at Los Alamos in collaboration with G. Kilcup, S. Sharpe and P. Tamayo. We end with 
a brief statement of the performance of our SIMD code on the CM5. 



1. LATTICE PARAMETERS 

All tlio results described in these two talks have 
been obtained using the following lattice param- 
eters. The 32^ gauge lattices were generated at 
/? = 6.0 using the combination 10 over-relaxed 
(OR) sweeps followed by 1 Metropolis sweep. We 
have stored lattices every 2500 OR sweeps. Quark 
propagators, using the simple Wilson action, have 
been calculated after doubling the lattices in 
the time direction, i.e. (32^ x 32 32^ x 64). 
This doubling was done in response to the initial 
hardware/operational constraints of the CM5 at 
LANL. Currently we arc generating 32^ x 64. Pe- 
riodic boundary conditions are used in all 4 direc- 
tions, both during lattice update and propagator 
calculation. Quark propagators have been calcu- 
lated using two kinds of extended sources - Wup- 
pertal and Wall - at k = 0.135 (C), 0.153 (S), 
0.155 (Ui). 0.1558 {U2), and 0.1563 {U3). These 
quark masses correspond to pseudoscalar mesons 
of mass 2800, 980, 700, 550 and 440 MeV re- 
spectively where we have used 1/a = 2.3GeV for 
the lattice scale. The three light quarks allow us 
to extrapolate the data to the physical u and d 
quarks, while the C and S k values are selected to 
be close to the physical charm and strange quark 
masses. The C (S) quark mass is slightly lower 
(higher) than the physical value. So far our sta- 
tistical sample consists of 28 configurations, each 
separated by 5000 OR sweeps. The goal is to an- 
alyze 100 configurations so the present analysis 
should be considered preliminary. Also, Sharpe 
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[4] and Bernard and Golterman [3] have pointed 
out that in the quenched approximation rj' loops 
give rise to unphysical terms in the chiral expan- 
sion for observables. These effects have not been 
included in the analysis when extrapolating quan- 
tities to the chiral limit. 

We find that the data for almost all oliservables 
show a significant curvature as a function of the 
quark mass. So we make both a linear and a 
quadratic fit in each case and quote final numbers 
from the fit with lower x^- More details about the 
fits are given in the individual sections. 

2. SPECTRUM 

We have analyzed three types of hadron cor- 
relators distinguished by the type of source/sink 
used to generate quark propagators. These are 
(i) wall source and point sink (WL), (ii) Wupper- 
tal source and point sink (SL), and Wuppertal 
source and sink (SS). With these three cases we 
are able to make a consistency check since the 
effective mass me//(t) converges to the asymp- 
totic value from below for WL correlators and 
from above for SL and SS correlators in all hadron 
channels. Fig. 1 illustrates, using the pion UiUi 
channel, noteworthy features of the convergence 
of m,pjf(f). In the case of WL correlators, the 
approach of meff{i) to its asymptotic value is 
monotonic but quite slow and one needs to fit 
for t > 15. The statistical errors in ■mr.ff(t) are 
smallest of the three cases. The SS correlators 
reach the asymptotic value at i « 8 and the time- 
slice to time-slice fiuctuations are of the order of 
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Figure 1. Comparison of the convergence of 
meff(t) for three types of pion (UiUi) correla- 
tors. SL data has been shifted by +0.05. 
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Figure 2. Quadratic extrapolation of A, nucleon 
and p mass data to the chiral limit. 
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the statistical errors. The SL correlators show 
correlated fluctuations in me//(t), wiggles with a 
period of about 12 time slices, which may be a 
feature of the Wuppertal source. Another pos- 
sibility, which we will check on the new set of 
32^ X 64 lattices, is that this feature is duo to 
doubling the lattice in the time direction. Inspite 
of these differences in me//(i), the bottom line 
is that all three kinds of correlators give consis- 
tent (within 1 (j) results for masses for all mesons 
and the error estimates are similar too. We give 
the average values for masses in lattice units in 
Table 1. 

To extrapolate the data to the chiral limit in 
the light quark wo fit as a function of m\ jrn^p and 
mq, where the latter is defined either as log(l + 
0.5(1/k— 1/kc)) or 0.5(1/k— 1/kc), or determined 
non-perturbativoly as discussed in Ref. [1]. Using 
a quadratic extrapolation of m\ data versus 
(as it fits all the points) we get 

= 0.15721(7). (1) 

This value is higher that given in Ref. [1], but 
consistent with QCDPAX result 0.15717(3). Sim- 
ilarly, using a quadratic fit, the extrapolated value 
of nip is (see Fig. 2) 

rap{ra^ = 0) = 0.34(4) (2) 



which gives a lattice scale of 
a-^{rap) = 2.26(27) GeV. 



(3) 



We shall, for simplicity, use l/a = 2.3 GeV. 

Wuppertal source quark propagators allow one 
to construct hadron correlators with non-zero mo- 
menta. In Fig. 3 wo show the pion propagator at 
4 different momenta. The signal in p = (1,0,0) 
and (1,1,0) correlators is quite good. Thus, as 
we show later, using 32"' lattices at /3 = 6.0 leads 
to a significantly improved signal in the calcula- 
tion of phenomenologically interesting quantities 
like semi-leptonic form factors and Bk- 

2.1. Baryons 

In our previous work[l] we had found that wall 
and Wuppertal sources give results that are sig- 
nificantly different. In Fig. 4 we show a typical 
example of the current status. We find that on 
the larger lattices WL and SL correlators give 
consistent results. The plateau in meff{i) is not 
very clean and any systematic difference between 
the two cases is an artifact of the fit range and/or 
the statistical errors. The numbers given in Ta- 
ble 1 for the case of degenerate quarks are the 
mean of SL and W L results. Extrapolating the 
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Figure 3. Effective mass plot for pion (UiUi) at 
different momenta. 
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data to tlie cliiral limit (Fig. 2) gives 

mjva = 0.44(20) 
mAfl = 0.59(40) . 



(4) 



however, before comparing these numbers to ex- 
perimental data the reader should note that we 
have not included the results of Labrenz and 
Sharpe[6] who predict additional non-analytic 
terms in the chiral behavior of quenched QCD. 

3. DECAY CONSTANTS 

Meson decay constants can be calculated using 
SL and SS correlators. We use the generic no- 
tation fps for pseudoscalar mesons and /V for 
vector meson. There are many different ways of 
combining these correlators to extract fps and 
/t/, some of which are described in Ref. [1]. We 
find that all methods give essentially identical re- 
sults, so we quote the mean value in Table 2. For 
the error estimate we give the largest individual 
error. We also find that the value of fps is in- 
dependent of the pion's momentum. For compar- 
ison wo give results using the local current with 
two different renormalization factors. Columns 2 
and 4 are using the Lepage-Mackenzie mean-field 
improved scheme[5] while columns 1 and 3 use the 



Table 1 

Masses in lattice units 



Nucleon 



cc 


1.219(2) 


1.232(02) 






cs 


0.857(2) 


0.884(03) 






CUi 


0.818(3) 


0.848(04) 






CU2 


0.803(4) 


0.833(05) 






CUs 


0.794(6) 


0.826(06) 






SS 


0.427(2) 


0.509(03) 


0.792(15) 0.854(19) 


SUi 


0.368(2) 


0.467(04) 






SU2 


0.344(2) 


0.451(05) 






SUs 


0.328(2) 


0.440(06) 








0.302(2) 


0.423(05) 


0.643(14) 0.721(22) 


U2U2 


0.240(2) 


0.389(09) 


0.572(12) 


0.668(27) 


U3U3 


0.193(2) 


0.364(12) 


0.526(16) 0.636(30) 


Table 2 








and in 


attice units 










-j^ y jnaive 


i//r'- 


CC 


0.126(3 


) 0.205(6) 


0.097(2) 


0.20(1) 


cs 


0.103(3 


) 0.137(4) 


0.121(3) 


0.20(1) 


CUi 


0.094(3 


) 0.122(4) 


0.116(3) 


0.19(1) 


CU2 


0.090(3 


) 0.116(4) 


0.115(4) 


0.19(1) 


CUs 


0.088(4 


) 0.113(5) 


0.114(4) 


0.18(1) 


SS 


0.090(3 


) 0.098(3) 


0.228(6) 


0.31(1) 


SUi 


0.083(3 


) 0.089(3) 


0.241(7) 


0.32(1) 


SU2 


0.080(3 


) 0.084(3) 


0.246(7) 


0.33(1) 


SU3 


0.078(3 


) 0.081(3) 


0.248(8) 


0.33(1) 




0.077(3 


) 0.080(3) 


0.261(7) 


0.34(1) 


U2U2 


0.071(3 


) 0.072(3) 


0.276(8) 


0.35(1) 


U3U3 


0.066(4 


) 0.067(4) 


0.285(9) 


0.36(1) 



naive v2k renormalization with Z,, = 0.57 and 
the 1-loop perturbative value for Za- In both 
cases we use a boosted value for the coupling, 
g'^ = 1.75. The results are displayed in Fig. 5. 

As shown in Fig. 5, data for f^^ is qualita- 
tively different for cases with degenerate versus 
non-degenerate quark masses. The large depen- 
dence on the renormalization constant (naive ver- 
sus improved) even for small nig suggests possibly 
large 0(«) corrections or effects of quenching in 
fp^- To get we can extrapolate the degener- 
ate combinations or all data with S and lighter 
quarks as they give qualitatively similar results. 
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Figure 4. Effective mass plot for nucleon with 
WL and SL correlators. 
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Using the mean-field improved normalizations we 
get (with all data) 



Ua = 0.052(5) 
Ua = 0.057(24) 



linear 
quadratic 



Both fits have comparable x^- From the linear 
fit wc got l/a = 2.5(3)Gey. For wc calculate 
the ratio fn/fK using S (Ui) as the strange quark 
and get the result 1.38 (1.52), where the change 
due to the light quark mass is « 1 in the last dec- 
imal place. Using ffc = l60MeV and the mean 
of the above two numbers (« physical s quark) 
we get fn = 232(20)Mey. 

4. Bk 

The major problem in the calculation of the 
kaon _B-parameter with Wilson fermions is the 
bad chiral behavior induced by the mixing of the 
AS = 2 4-fermion operator with wrong chirality 
operators. This mixing arises as a result of the 
r-term in the action and has been calculated to 
1-loop in perturbation theory. In Ref. [2] it was 
shown that the lattice artifacts completely over- 
whelm the signal when using the 1-loop improved 
operator, however, by calculating the matrix el- 
ement at different values of momentum transfer 



Figure 5. Meson decay constants fps and f^l 
as a function of m^//?i^. Data for fy^ has been 
shifted by +0.25 for clarity. 
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one can remove most of these. The main limita- 
tion of that calculation was that the spatial lattice 
size was small (16^), consequently the momen- 
tum gap (27r/16) was large, and the signal in the 
non-zero momentum correlators was marginal. 
Both these limitations have been addressed on 
the present set of lattices and as a result the sta- 
tistical quality of the data is markedly improved. 

In Fig. 6 we show the behavior of the lattice 
Bk, evaluated at momentum transfer p = 27r/32, 
as a function of the time slice at which the op- 
erator is inserted. The data shown are for two 
different kaon sources, i.e. tps^si'd find V'.s'7475'0(i 
as they converge from opposite direction. The 
statistical errors are small and the two fits agree 
within errors. The results for p = 27r/32 x (1, 1, 0) 
and p = 27r/32 x (1,1,1) are also reliable. Full 
analysis will be presented elsewhere Ref. [7]. 

In Table 3 we present preliminary results for 
the _B-parameters using perturbatively improved 
operators evaluated with a boosted = 1.75 
and //,« = 1.0. In column 1 (2) we give the raw 
number for zero (one) unit of momentum trans- 
fer. The final values of Bk after removing the 
leading bad chiral contributions by the method 
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Figure 6. Bk at momentum transfer p = (1,0,0) 
as a function of t at which the operator is inserted. 
The two data sets correspond to the two kinds of 
operators used to produce the kaons. 
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of momentum subtraction (Ref. [2]) are given in 
column 3. In columns 4 and 5 we give the re- 
sults for the two left-right operators O7 and Og 
(defined in Eof. [2]). 

The data in columns 1 and 2 show that the 
statistical errors have been reduced to the level of 
a few percent. The reason that the errors in the 
final _B-parameters are large for heavy quarks is 
mainly because E(p) « in, while for light quarks 
the statistical fluctuations are largo. 

The numbers extrapolated to the kaon mass are 
given in the last two rows in Table 3. We use two 
methods for this: in the first we keep S mass fixed 
and extrapolate in the light quark mass using 
S, Ui, U2, U3. In the second, the two quarks are 
taken to be degenerate. Both methods give sim- 
ilar results for all three operators. This suggests 
that the SU(3) breaking does not give a large con- 
tribution to these _B-parameters. Also, there is 
little evidence of unphysical contributions from 
rj' loops which are expected to be significant near 
the chiral limit in the quenched approximation[8]. 

The values of Bk, B-j and Bg are significantly 
lower than those found in Ref. [2]. We believe 
that this is due to using a larger lattice and 



Table 3 
B^. D-j and 







Bk 




B7 


Bs 




(p = 0) 


(p= 1) 


subtr. 






cc 


0.95(3) 


0.95(4) 


1.32(37) 


0.97(3) 


0.97(3) 


cs 


0.83(2) 


0.84(3) 


1.09(21) 


0.94(3) 


0.95(3) 


CUi 


0.79(2) 


0.80(2) 


0.94(19) 


0.92(3) 


0.94(3) 


CU2 


0.77(2) 


0.77(2) 


0.76(20) 


0.91(3) 


0.93(3) 


CLh 


0.76(2) 


0.75(3) 


0.50(24) 


0.89(3) 


0.92(3) 


SS 


0.50(1) 


0.53(1) 


0.82( 7) 


0.87(2) 


0.91(2) 




0.36(1) 


0.41(2) 


0.72( 6) 


0.84(1) 


0.89(1) 


SU2 


0.29(1) 


0.34(2) 


0.66( 7) 


0.83(1) 


0.88(1) 


SU3 


0.23(1) 


0.28(2) 


0.60( 8) 


0.82(1) 


0.87(1) 




0.13(1) 


0.21(2) 


0.62( 6) 


0.81(1) 


0.86(1) 


U2U2 


-.26(2) 


-.08(3) 


0.48( 8) 


0.76(1) 


0.82(1) 




-.78(3) 


-.42(4) 


0.37(11) 


0.72(1) 


0.78(1) 


Ext 1 






0.48 


0.80 


0.85 


Ext 2 






0.43 


0.74 


0.80 



a smaller p. Very accurate and reliable results 
for Bk, obtained using staggered fermions, have 
been presented at this conference by S. Sliarpe[8]. 
Wilson results are consistent with these but have 
larger systematic errors as a result of mixing with 
operators having the wrong chiral behavior. By 
improving the quality of Wilson fermion results 
one can, by comparison with staggered results, 
check our understanding of lattice artifacts. 

5. SEMILEPTONIC MESON FORM FACTORS 

We have calculated the semileptonic pseudo- 
scalar form factors of the D meson, both for 
D — ^ Kev and D — )• ttcp. In calculating these 
form factors we have held the D-meson at zero 
spatial momentum, and varied the momentum of 
the other meson from to 77/8. This provides a 
large enough range (both positive and negative) 
in the invariant mass of the leptonic subsystem 
(—Q^) to test the pole dominance hypotheses. To 
control systematic effects, we have done the calcu- 
lation using three different transcriptions (called 
'local", 'extended' and 'conserved") of the current. 
The details of the method are given in Ref. [9] and 
a complete analysis for the present set of lattices 
will be presented in Ref. [10]. 
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Tabic 4 

Tlio meson form factors extrapolated to = 

f^jo) fm /o^(Q) /o"(o) 

Ui 0.833(21) 0.787(24) 0.808(14) 0.735(14) 
U2 0.837(24) 0.761(32) 0.808(15) 0.699(16) 
Us 0.834(27) 0.737(41) 0.806(17) 0.670(18) 



The renormalization of all three currents is 
carried out using the Lepage-Mackenzie[5] mean 
field improved perturbation theory. We find that 
the values obtained using the three currents are 
very close. (See Fig. 7 for an example of con- 
sistency between different currents.) We notice 
that the statistical errors are very small when 
the hadrons carry zero or one unit of momentum, 
and some systematic difference between the local 
and the non-local currents become visible. For 
p = (1, 1, 1) and (2, 0, 0) there is no clear plateau 
in the data and the analysis is very sensitive to 
the fit range. 

Except when p = (1,1,1), the data are con- 
sistent (see Fig. 7 for an example) with the 
hypothesis that the form factor is dominated by 
the nearest pole with the right quantum numbers 
(D* and D* respectively for and /^). Assum- 
ing pole dominance we can extract f(Q^ = 0), 
and the results are summarized in Table 4. To 
within statistical errors the numbers are indepen- 
dent of the light quark mass for and there is 
a slight decrease in /'^. We expect « 10% change 
in /"^ on extrapolating the data to the chiral 
limit and to the physical values of charm quark 
mass. Our final numbers for f:^ are close to the 
phenomenological value, « /(f and f:^ > f^. 



Figure 7. Comparison of I//0 obtained using the 
three vector currents. The data are for the case 
where the spectator quark is C/3. The fits show 
comparison with pole dominance using the mea- 
sured lattice value for the scalar mass, to^^^. 
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35 megaflops/node. These timings include all I/O 
and setup overhead. For the present set of physics 
objectives we need 18000 node-hours to generate 
and process one configuration. 

We gratefully acknowledge the tremendous 
support provided by the ACL at Los Alamos. 
These calculations have been done as part of the 
DOE HPCC Grand Challenge program. 
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